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We report on the early time ordering in a nematic liquid crystal subjected to a sudden change in 
external ac electric field. We compare time evolution for two different initial states of electrocon- 
vection. Electroconvection is a highly driven state of a nematic liquid crystal involving convect ive 
motion of the fluid and periodic variations of the molecular alignment. By suddenly changing either 
the voltage or the frequency of the applied ac field, the system is brought to the same thermody- 
namic conditions. The time ordering of the system is characterized by the evolution of features of the 
power spectrum, including the average wavenumber, total power, and shape of the the power spec- 
trum. The differences between the two classes of quenches are discussed, as well as the possibility 
of scaling behavior during this initial phase of domain growth. 
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Understanding systems driven far from equilibrium re- 
mains one of the outstanding challenges of contemporary 
physics. At the heart of the issue is the lack of any sin- 
gle principle that is equivalent to the minimization of 
free energy that is applicable in thermodynamic equilib- 
rium A subset of this larger issue is the question of 
the transition between states of a system after a sudden 
change in an external parameter, or a quench. The gen- 
eral question of the behavior of systems after a quench is 
often referred to as phase ordering or coarsening |2|, as 
domains of the new steady state of the system order and 
grow. 

Classically, phase ordering has been studied in the con- 
text of the transition between two equilibrium states. In 
this case, minimization of a free energy plays a central 
role in understanding the ordering process |2|. More re- 
cently, transitions between driven states of a system have 
gained interest 0, H H H 0, HH ED, El El In this case, 
one considers the ordering of steady-states of a driven 
system after the driving force has been changed. The free 
energy does not play a role. However, there is growing 
evidence that the dynamics of topological defects govern 
the phase ordering in both equilibrium and nonequilib- 
rium systems 

H El E3 One system that is useful for 
studying quenches in driven systems is electroconvection 

mum 

Electroconvection occurs when a nematic liquid crystal 
is placed between two plates and an ac electric voltage is 
applied perpendicular to the plates 0, ^| . A nematic 
liquid crystal consists of long, rod-like molecules that on 
average are aligned along a fixed axis, which is referred to 
as the director. In the absence of electroconvection, the 
director is spatially uniform. Above a critical voltage, the 
director develops a periodic spatial variation, and there is 
an associated periodic flow, or convection rolls. The most 
common geometry for electroconvection is to have the di- 
rector parallel to the glass plates initially. An interesting 
case exists when the director is initially perpendicular to 



the plates, also known as homeotropic alignment [TfjlT^ . 

An important parameter that is used to characterize 
nematic liquid crystals is the dielectric anisotropy [TEj . 
For an anisotropic material, such as a nematic liquid 
crystal, the dielectric constant is a tensor. The dielec- 
tric anisotropy is the difference between the dielectric 
constant for the case where the electric field and direc- 
tor are perpendicular and the case where they are par- 
allel. For materials with a negative dielectric anisotropy 
in a homeotropic configuration, the initial transition as a 
function of applied voltage is the Freedericks transition, 
which is followed by the transition to electroconvection 
[Trl ll^ . The Freedericks transition corresponds to the 
director developing a tilt relative to its initial alignment 
perpendicular to the plates. This state is an equilibrium 
state in which the angle of tilt is fixed for a given value 
of the voltage. However, there is a degeneracy due to the 
fact that the director tilt can adopt any azimuthal an- 
gle. Regions with the same azimuthal angle are usually 
referred to as Freedericks domains. Once the Freedericks 
transition has occurred, there is a component of the di- 
rector parallel to the plates. Therefore, there is a second 
critical voltage at which electroconvection will occur. 

The Freedericks transition is independent of frequency; 
whereas, the critical voltage for electroconvection is fre- 
quency dependent. This allows for exploration of two dif- 
ferent types of quenches to the same equilibrium state. 
With different starting points, either a rapid change 
in frequency or voltage can bring the system to the 
same final frequency and voltage, and hence the same 
Freedericks state. However, the system evolution to reach 
this final state is different. For a quench down in volt- 
age, the average director tilt angle has to relax to the 
correct value. For a quench in frequency, the average 
tilt angle of director should already have the appropriate 
value. In both cases, the pattern present in the director 
field, charge distribution, and flow field, due to electro- 
convection have to relax. 
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In this paper, we report on a comparison of these two 
quenches. We focus on the initial dynamics from the time 
of the quench until the point where the spatial patterns 
are similar. The system's evolution is characterized by 
a number of measures based on the power spectrum of 
the images. It should be noted that it is also interest- 
ing to consider the late-time evolution of the system. As 
the Freedericks state is an equilibrium state, the late- 
time evolution should be the same for both systems Q. 
However, large-scale imperfections in the system compli- 
cate studies of the late-time dynamics, for which large 
aspects ratios are critical. Therefore, the question of the 
late-time behavior will be the subject of future work. 

For these experiments, homeotropic cells of 
the nematic liquid crystal N4 were used. The 
liquid crystal N4 is a eutectic mixture of the 
two isomers of 4-methoxy-4'-n-butylazoxybenzene 
(CH 3 - C 6 H 4 - NON - C 6 H 4 - C4H9 and 
CH3OC6H4-NNO-C6H4-C4H9). It was ob- 
tained from EM Industries (a Merck company), now 
EMD Chemicals Inc. [l^ The N4 was used with- 
out further purification. The method for obtaining 
homeotropic alignment is described in detail in Ref. [HI . 
Briefly, a surfactant coating is made on ITO coated 
glass using a Langmuir-Blodgett technique. A mixture 
of 43% of N4 and 57% arachidic acid (C20) diluted in 
chloroform is spread at the air- water interface, forming 
a monolayer that was compressed to 10 mN/m. We used 
C20 obtained from Sigma- Aldrich with a quoted purity 
of > 99%. It was used without further purification. The 
pressure was held constant while ITO glass is coated 
with 10 layers of the Langmuir monolayer. After the 
coating of the surface, the glass is baked in an oven at a 
temperature of 50 °C. A 25 jam mylar spacer is placed 
between the two ITO glasses. Two opposite sides of the 
cell are sealed with epoxy. The cell is filled with N4 by 
capillary action. The two remaining sides of the cell are 
sealed with 5 minute epoxy. 

The details of the experimental apparatus are de- 
scribed in Ref. [2(J. Essentially, there is a temperature 
control stage that holds the sample at a constant temper- 
ature within ±5 mK. The sample is illuminated from be- 
low and imaged from above using standard shadowgraph 
techniques. A pair of crossed-polarizers (one below the 
sample and one above the sample) allow for simultaneous 
imaging of electroconvection patterns and of Freedericksz 
domains. The use of ITO coated glass allows for the ap- 
plication of an ac voltage to the sample. 

The general state diagram for homeotropic samples 
of N4 is reported in Ref. 0. There is the expected 
Freedericks transition at approximately V c f = 9.3 V, in- 
dependent of applied frequency. The transition to elec- 
troconvection is frequency dependent. The solid black 
horizontal and vertical lines with triangular endpoints in 
Fig. 1 illustrate the two types of quenches that are re- 
ported on in this paper. The solid curve is the transition 
to electroconvection, and the dashed black curve is the 
critical voltage for the Freedericks transition (both curves 
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FIG. 1: A portion of the phase diagram for homeotropic N4 as 
reported in Ref. [19| . The dashed curve is the critical voltage 
for the Freedericksz transition. The solid curve is the critical 
voltage for the transition to electroconvection. The horizontal 
and vertical lines connect the starting and ending point of the 
two quenches. 

were reported in Ref. ) . Both quenches are selected to 
have the same final equilibrium conditions of applied volt- 
age and frequency within the regime where Freedericks 
domains exist, but electroconvection does not. The two 
starting points represent different initial states for the 
system. The quench labelled A starts from a chaotic elec- 
troconvection state and uses a change in voltage. Phys- 
ically, this means the tilt of the director must change as 
part of the quench. The path labelled B starts from a reg- 
ular state of electroconvection and involves a change in 
frequency. In this case, the distance from the Freedericks 
critical voltage is not changed. Therefore, the details of 
the initial evolution should be different. 

For discussing these quenches, it is important to de- 
fine a few parameters. We will use e = (V/V c ) 2 — 1 to 
refer to the distance from the electroconvection critical 
voltage (V c ) at a fixed frequency. Therefore, since all of 
the quenches end at a frequency / = 5000 Hz and a rms 
voltage of V = 13.0 V, the final point for the quenches 
is e = —0.156 and / = 5000 Hz. The two starting points 
are e = 0.618 with / = 5000 Hz and e = 0.490 with 
/ = 25 Hz. 

To provide a framework for the two quenches, we have 
characterized the final state both by making quenches 
from below the Freedericks critical voltage and by slow- 
ing stepping to the final state. These studies involving 
increasing the voltage established the existence of large- 
scale spatial inhomogeneities in the Freedericks domains. 
Essentially the same spatial pattern is obtained either 
by small steps in voltage or by considering the late-time 
state after a large change in voltage. The spatial inho- 
mogeneities are most likely caused by small defects in 
the aligning layer that result in the pinning of defects 
in the director organization. Because of the large-scale 
nature of the inhomogeneities, we do not expect them 
to play a significant role in the early-time evolution of 
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the system. However, they do prevent detailed studies 
of the late-time ordering at this point, and the influence 
of these domains can not be completely ruled out for the 
early time evolution. 

For each of the quenches, we measure a series of 128 
images taken 1 s apart. We characterized the images us- 
ing the spatial Fourier transform of each image. We focus 
on the square of the modulus of the Fourier transform of 
each image (the power spectrum, S(k,6), where k is the 
wavenumber of interest). The power spectra of 20 jumps 
are averaged image by image to improve the statistics. 
There is a 15 minute wait between jumps. Using 5(fc, 0), 
we measure the total power 

p = Jo Jo°° S(k,6)kdkd0 = 
2n J °° S(k)kdk, where S(k) is the azimuthal average of 
S(k,9). From this measure, we confirmed that the total 
power is essentially constant throughout the time evolu- 
tion, allowing comparison between features of the power 
spectra at different times. The average wavenumber, 



< k >- 



27r/ °° S(k)k 2 dk 



(1) 



is used as a measure of the typical length scale in the sys- 
tem. The rate of growth of this length scale can be com- 
pared for the two quenches using the time-dependence of 
the average wavenumber. Because we are looking at early 
time dynamics, it is not necessarily expected that a scal- 
ing regime will exist. However, we tested for scaling of 
the power spectrum of the form S(k, t) = L d g(kL\ where 
L represents a typical length scale for the system and d is 
the spatial dimension of the pattern. For our case, d = 2, 
and we take L =< k > _1 . Finally, in the study of phase 
ordering, Porod's Law @, E[ states that in the limit of 
large wavelength (kL » 1), S(k,t) ~ L~ n k~^ d+n) . In 
this case, d = 2 and n is the dimension of the vector 
order parameters that describes the system. For a ne- 
matic liquid crystal, one expects n = 2. In this sense, 
Porod's Law relates the behavior of the long wavelength 
tail of the power spectrum to the dominate defects in the 
system Q. 




FIG. 2: Series of images for the frequency (a - d) and voltage 
(e - h) quenches. The white bar in the images represents 
1 mm. The images are taken at 1 s, 4 s, 32 s, and 128 s after 
each quench. 

Figure 2 presents the typical time evolution of the sys- 
tem after a quench in voltage. Images 2(a) - (d) are 
for a frequency quench, and images 2(e) - (h) are for a 



voltage quench. The images are taken at 1 s, 4 s, 32 s, 
and 128 s after the quench and illustrate the ordering 
that occurs. The images suggest that the voltage quench 
evolves slower than the frequency quench. To quantify 
this difference, we used the time dependence of < k >. 
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FIG. 3: Plot of the average wavenumber versus time for each 
of the quenches. The frequency quench is represented by 
squares and the voltage quench is represented by triangles. 
The solid line in each case is a fit to a power law. The ex- 
ponent for the frequency quench is -0.47 and for the voltage 
quench, it is -0.37. 

The main result of this paper is shown in Fig. 3. Plot- 
ted in Fig. 3 is the average wavenumber, /c, versus time 
for both quenches. The most dramatic feature of this 
result for the frequency quench is the agreement with 
power law growth over the entire time range. The volt- 
age quench is also consistent with power law growth for 
the time range plotted in Fig. 3. By comparing the ex- 
ponents in the power law (-0.47 for the frequency quench 
and -0.37 for the voltage quench), one finds that the sys- 
tem orders faster after a frequency quench than it does 
after a voltage quench. 

The existence of power law behavior for the evolu- 
tion of the average wavenumber suggests that the sys- 
tem might be in a scaling regime. The limited time of 
observation makes it difficult to establish this with any 
certainty. However, as a test of scaling, we considered the 
behavior of the full power spectra. The results for scaling 
the azimuthally averaged power spectra are shown in Fig. 
4. The results are consistent with a scaling of the power 
spectra, but are not conclusive. Perhaps more important 
than the overall behavior of the spectra is the behavior 
at large wavenumber. This is shown in the insert of Fig. 
5 for a time 16 s after the quench. The long wavelength 
part of the power spectra is consistent with a power law 
that is independent of the quench type: S(k) ~ k~ 4 for 
large k. This result is in agreement with the expected 
value for ordering of a nematic liquid crystal in two di- 
mensions 0- 

To summarize, the two quenches started from different 
values of the average tilt of the director. Because the 
average wavenumber was consistent with a power law, 
the different rates of phase ordering could be quantified. 
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FIG. 4: Plot of the scaled power spectra versus the average 
wavenumber for the frequency quench (solid symbols) and for 
the voltage quench (open symbols). The power spectra for 2 s, 
32 s, 64 s, and 128 s after the quench is shown for each case. 
The insert shows the long time tail for 16 s after the frequency 
quench (squares) and the voltage quench (triangles). The 
solid line represents k~ A and is provided as a guide to the 
eye. 

The frequency quench (with the same initial and final 
tilt) exhibited faster ordering than the voltage quench. 
This is consistent with the different intrinsic time scales 
for electroconvection. The two main relaxation times are 
the charge relaxation time (on the order of 10 -3 s for this 
system) and the director relaxation time (on the order of 
1 s for this system). Because the voltage quench should 
be dominated by director relaxation at early times as the 



system equilibrates to the equilibrium tilt angle, it is not 
too surprising that this quench has slower dynamics. 

Of interest for future work is the behavior of the ex- 
ponents that correspond to the power law growth of the 
average wavenumber. For this system, the order param- 
eter is not conserved, so one would expect an exponent 
of 1/2. This work was a focused study that compared 
frequency quenches and voltage quenches. For the fre- 
quency quench, the scaling of the average wavenumber 
was close to 1/2 (0.47), but for the voltage quench, it 
was significantly smaller (0.37). It should be noted that 
for the voltage quench, by considering different time win- 
dows (0.1 s, 0.5 s, and 1 s), we observed apparent expo- 
nents for the evolution of < k > of 0.28, 0.35, and 0.37, 
respectively. For the frequency quenches, the exponents 
were essentially constant. This suggests two things. The 
time dependence of the evolution of the two systems may 
become similar in the long-time limit, if the trend for 
the voltage quench continues. If that turns out to be 
true, both systems obey the expected growth laws for an 
equilbrium system at late times. However, more uniform 
systems and longer times are needed to confirm this be- 
havior, with a focus on the voltage quenches for which 
the behavior is the most dramatic. Additionally, the im- 
pact of the type of quench can be further explored using 
mixed quenches where both the frequency and voltage is 
varied. 
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